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This paper describes the calculation method for unsteady state conditions in the secondary air systems in gas tur- 
bines. The 1D-3D-Structure coupled method was applied. A 1D code was used to model the standard components 
that have typical geometric characteristics. Their flow and heat transfer were described by empirical correlations 
based on experimental data or CFD calculations. A 3D code was used to model the non-standard components that 
cannot be described by typical geometric languages, while a finite element analysis was carried out to compute 
the structural deformation and heat conduction at certain important positions. These codes were coupled through 
their interfaces. Thus, the changes in heat transfer and structure and their interactions caused by exterior distur- 
bances can be reflected. The results of the coupling method in an unsteady state showed an apparent deviation 
from the existing data, while the results in the steady state were highly consistent with the existing data. The dif- 
ference in the results in the unsteady state was caused primarily by structural deformation that cannot be predicted 


by the 1D method. Thus, in order to obtain the unsteady state performance of a secondary air system more accu- 


rately and efficiently, the 1D-3D-Structure coupled method should be used. 
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Introduction 


The earliest exploration of one-dimensional (1D) fluid 
network calculations dates to 1936, when Cross[1] pro- 
posed the application of network iterations for the calcu- 
lation of steady flow. In later decades, the theory and 
methods pertaining to this subject were integrated with 
experiments and improved gradually. A number of insti- 
tutes and companies attempted to associate the theory of 
a flow network with engineering and developed a useful 
tool for complicated fluid networks. Yet it was not until 
the 1990s that the method of the flow network was truly 
adopted in the secondary air system of gas turbines. Kutz 
and Speer [2] reviewed a computer program for the si- 
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mulation of a secondary air system. An algorithm for the 
solutions to the resulting nonlinear governing equations 
was introduced, and the loss correlations of the typical 
flow elements were also discussed. The results of the 
program were verified against test data with strong 
agreement. Muller [3] used the Newton-Raphson method 
to calculate the nonlinear equations for the set of un- 
known parameters, including total pressure, total temper- 
ature, mass flow rate, and others. Furthermore, the net- 
work was coupled within a thermo-mechanical finite 
element model of an engine’s secondary air system 
structure with the free finite element software CalculiX. 
The results of the coupled investigation showed good 
agreement with test data, but heat transfer data were still 
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derived on the basis of empirical correlations. Past re- 
search on secondary air systems focused mainly on 
steady state conditions [4]. However, the gas turbines in 
airplanes always work in unsteady state conditions that 
are caused primarily by the following: 

1. Starting and stopping of the engines. 

2. Change in rotational speed and power output. 

3. Change in the altitude of the engine. 

4. Fluctuations in the pressure and temperature at the 
outlet of the compressor. 

5. Storage or dissipation of mass and energy in the 
system. 

Moreover, the deformation of hot components may 
occur in concert with these factors and influence the flow 
and heat transfer inside the system. These phenomena 
can result in the ingress of gas into the cavities, which 
can destroy the pressure balance in the system. Hot parts 
may become overheated, and the efficiency of sealing 
can decrease as a result. Thus, whether or not the sec- 
ondary air system can supply sufficient air cooling with 
proper pressure to prevent the system from failing be- 
comes extremely important. In other industrial fields, 
such as petroleum, natural gas, internal-combustion en- 
gines, and pulse detonation engines, among others, re- 
search on unsteady state performance, such as water 
hammer, gas hammer, fuel injection, dynamic drainage, 
and pulse explosion, was carried out extensively by in- 
troducing the concept of “pressure wave and disturbance” 
into the calculations. Furthermore, methods pertaining to 
characteristics, finite difference, finite volume, and finite 
element were developed. 

Since the 1960s, Benson and his team have conducted 
research on unsteady state conditions in internal combus- 
tion engines and the calculation of characteristic sets of 
equations. In the 1970s, his team finished the UMIST 
program and published a monograph [7] that summarized 
the research on the method of characteristics. For a long 
period, the method of characteristics was the mainstream 
technique in 1D research of unsteady state conditions in 
internal combustion engines and in other industrial fields. 
With the development of computational techniques, the 
methods of finite difference, finite volume, and finite 
element have been applied more often, but the develop- 
ment of boundary equations is still based primarily on the 
theory of characteristics [8]. The traditional 1D unsteady 
state methods can predict the unsteady effectiveness of 
the accumulation or dissipation of mass and heat inside 
the system, but for the structural deformation and result- 
ing changes in flow, the 1D-Structure coupled method 
must be applied. For certain special and complicated 
structures that are difficult to describe with 1D correla- 
tion, the 3D tool is used and finally coupled together with 
1D and finite element analytical tools. 
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Coupling Mechanism of Multi-physical Fields 


For example, consider a two-fold physical phenome- 
non. When two different systems have their own govern- 
ing equations, coupling is realized by the interaction of c 
and d. 


f (a, b, c, d)=0 
and 
g (c, d, e, N=0 


Code A 


f (a, b, c, d)=0 


Coupling Zone 


Fig. 1 Coupling of two different codes 


The coupling zone indicates different physical boun- 
daries for different codes. For the 1D and 3D flow fields, 
the zone indicates the interface between them; for the 
flow field and solid body, the zone denotes a real wall 
surface. Coupling methods include direct, sequential, and 
synchronous coupling [9]. The method of direct coupling 
contains only one analysis and uses a coupling unit with 
multi-field free degrees. Coupling is carried out by cal- 
culating the cell matrix or load vector that contains the 
relative parameters. The method of sequential coupling 
contains more than two sequential analyses, and each 
pertains to one type of physical field; coupling is carried 
out by interacting the results of one analysis with another. 
The method of synchronous coupling indicates that dif- 
ferent fields are calculated in coincidence with each other, 
thus satisfying the boundary equations of each. This is 
strong coupling in the physical sense. 


1D-3D Coupling 

For the secondary air system in an aero-engine, an 
absolute 1D analysis has the following difficulties in pre- 
dicting the performance of the system: 

1. Some non-standard flow structures cannot be de- 
scribed with existing mathematical models, or the error is 
extremely large. 

2. Performance data of different components are un- 
certain. Although they can be derived from experiments 
or CFD simulations, the real-time effect over a transient 
course cannot be reflected. 

The joint simulation of 1D and 3D fluid dynamics 
shows the data transfer and iteration at the interface be- 
tween the 1D and 3D flow networks. In a secondary air 
system, the accuracy of cavity modeling is rather impor- 
tant, while the structure of the cavity is always compli- 
cated and variable, and thus 1D empirical correlations 
cannot describe the flow and course of heat transfer pre- 
cisely or are limited to certain confined working condi- 
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tions. Therefore, 1D-3D coupling should be carried out to 
obtain more reliable results. A 3D flow dynamics tool is 
used to model complicated cavities or other components, 
while normal parts or discrete losses are modeled with 
the 1D tool. Heat and flow boundary conditions are 
transferred between these two types of tools for flow 
dynamics. In this way, the unsteady state performance of 
a secondary air system can be derived. Figure 2 shows 
the scheme of the data transfer between Flowmaster, a 
commercial 1D flow dynamics code, and CFD code in 
which the interface is specified artificially. Mass flow 
rate, total pressure, total temperature, etc., are transferred 
through the two sides of the interface, and the results of 
the last iteration are transferred to the next iteration as 
boundary conditions [10]. 

The methods of data transfer at the interface are listed 
in Table 1. 


eee Flow direction 
ois » Data from Flowmaster to CFD 
desea > Data from CFD to Flowmaster 


j Three-dimensional F] = 
J ` CFD Model 7 


Flowmaster be: CFD be: CFD be: Flowmaster be: 
Pressure source Mass flow inlet Pressure outlet Flow source 


Mass Flow Rate ----- > 
---- Total Pressure 


Mass Flux Rate -----» 
---- Total Pressure 


Fig. 2 Data transfer between Flowmaster and CFD 


Table 1 1D-3D data transfer on the interface 
1D to 3D 3D to 1D 
Extensive Quantities wags A p=YV, 
(flowrate, velocity, etc.) A i= 
Intensive Quantities v,=p p= < v A 
(pressure, temperature, etc.) io A 


In Table 1, A represents the boundary surface area of 
the 3D model, and A; represents the unit area of each cell 
mesh on the boundary of the 3D model. Also in the table, 
p represents the boundary value of the 1D model, and v; 
represents the unit value of each cell mesh on the boun- 
dary surface of the 3D model. For data transfers from the 
1D to the 3D model, the mass flow rate is divided ac- 
cording to the percentage area of each cell face, while the 
pressure and temperature remain the same on each cell 
node. 


Fluid Solid Coupling 

In essence, all of the unknown variables in the fluid 
field and solid field should be calculated at the same time 
within one timestep for fluid solid coupling, or so-called 
strong coupling, although weak coupling is another me- 
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thod of fluid solid coupling. Zhang and Hisada of Tokyo 
University [11] distinguished strong coupling and weak 
coupling theoretically, compared their calculating effi- 
ciency, and then, on this basis, used the single physical 
equation as a starting point and clarified the essence of 
strong coupling, weak coupling, and one-way coupling. 
Establishment and Condensation of Coupling Equa- 
tions 


When solving the problems of fluid solid coupling by 
the finite element method (i.e., dispersing the fluid and 
solid fields by the finite element method), the boundary 
conditions are realized during the combination of ele- 
ments. If AX is the variable increment in the coupling 
system, according to the difference in the physical field 
where the nodes lie, then the simplified coupling equa- 
tions can be obtained, as shown below. 


Af Af- 0 |(axf RÍ 
Adi Adc +Aéc Ady || AX E RË (1) 
0 Aic Ar || AX; R; 


Where the superscripts f and s represent the fluid field 
and solid field, respectively, and the subscripts C and I 
represent variables on the coupling interface and interior 
nodes, respectively. AS and Af are the equivalent mass 


matrices of the coupling system. AX i ; AXE, and 


AX; are the unknown number vectors on the nodes in 


the fluid field, coupling interface, and solid field, respec- 
tively. RÍ, RË, and R are the external force vectors in the 
fluid field, coupling the interface, and solid field. There 
are two ways to solve the static force condensation in Eq. 
(1): first, eliminate the interior variables inside the fluid 
field and then calculate the solid field in combination 
with the interface, or revise the interior variable vectors 
inside the fluid field; second, eliminate the interior va- 
riables inside the solid field, and then calculate the fluid 
field in combination with the interface or revise the inte- 
rior variable vectors inside the solid field. 


Solution of Coupling Equations 


Based on the uniform coupling equations and the me- 
thod of condensation, strong coupling equations, coupl- 
ing variables in linear coupling equations, and the revi- 
sion of the variable vectors can be solved. Strong coupl- 
ing must deal with transferred information. According to 
the difference in boundary information between the fluid 
and solid fields, the interface can be classified as the es- 
sential boundary and natural boundary. The essential 
boundary defines the kinematic conditions, such as ve- 
locity, displacement and other boundary information, 
while the natural boundary defines the external force 
conditions, such as pressure, pressure force, and so on. 
Unlike strong coupling, structure and fluid solvers are 
carried out alternately in weak coupling. According to the 
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condensation method of coupling equations, there are 
four types of weak coupling methods: 

1. The common solution of structure and interface 
(WCS): Take the interface as the natural boundary for the 
solid field calculation and the essential boundary for the 
fluid field calculation. 

2. The common solution of fluid and interface (WCF): 
Because the free degree of the fluid field is much larger 
than that of the solid field, in order to solve the weak 
boundary conditions by WCS and additional huge time 


consumption of the mass matrix Afg calculation, the 


method of WCF is adopted. The interface is taken as the 
natural boundary for the fluid field calculation and as the 
essential boundary for the solid field calculation. 

3. The common solution of fluid, solid, and interface 
(WCFS): The interface is taken as the natural boundary 
for both fluid and solid fields. 

4. One-directional coupling: In some cases, the de- 
formation of the solid is rather slight, so the influence on 
the fluid can be ignored. This method can be taken as the 


simplified WCS method when the mass matrix Alia =0. 
1D-3D-Structure Coupling through MpCCI 


For secondary air systems in gas turbines, there are 
always complicated flow structures, such as cavities with 
multi-connections. Because the flow and heat transfer in 
the cavities influence the system to a great extent, the 
model of the cavities must be very accurate. A 1D model 
can model the cavity efficiently and simply, but the ac- 
curacy is mostly due to the embedded empirical correla- 
tions; thus, there exists much more uncertainty when this 
is applied to even more complicated structures. In order 
to obtain more reliable results for secondary air systems, 
this paper introduces a new method of coupling that inte- 
grates Flowmaster, Fluent and Abaqus through MpCCI. 
Fluent is used to model such complicated components as 
cavities, Flowmaster is used to model the standard struc- 
tures and discrete losses in the system, while Abaqus is 
imported to realize the bidirectional coupling of the flow 
network and structural network. As a result, the influence 
on the flow of sealing clearance changes caused by 
structural deformation, together with the influence of 
flow and heat transfer on the structure, can be taken into 
account. Furthermore, the method of 1D-3D-Structure 
coupling can eliminate the influence on the results of 
human factors when defining the boundary conditions at 
the interface of the fluid and solid fields artificially. 


1D-3D-Structure Coupling of a Secondary Air 
System 


Figure 5 shows the scheme of the flow network of the 
secondary air system in one type of gas turbine. Com- 
pressed air is imported into the conical stationary cavity 
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through a circular array of foil nozzles and orifices. Next, 
the air flows into an annulus passage with an inner wall 
that rotates through two circular arrays of orifices. Fol- 
lowing that, the air is divided into three branches. One 
branch discharges through a straight-through labyrinth 
seal at the left end of the passage. A second branch flows 
into the rotor-stator structure in front of the rotating disk 
through a straight-through labyrinth seal at the right end 
of the passage. Because of the centrifugal effect of rota- 
tion, this portion of the air is pumped out of the cavity 
through the rim seal and mixed with the mainstream gas. 
The third branch flows into the central rotating passage 
through a circular array of orifices, and then flows suc- 
cessively across round-sectional long orifices, discrete 
losses, and the rotor-stator structure at the back of the 
rotating disk, where it finally joins the mainstream gas 
through the rim seal structure. 

The whole system contains standard structures, such 
as rotating orifices, straight-through labyrinth seals, and 


CFD Solver Finite Element Solver 
K= F,(t) MOCK AD: 
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Fig. 3 Scheme of Fluid Solid Coupling 
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Fig. 5 Scheme of the flow network of the secondary air sys- 
tem in one type of gas turbine 
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annulus passages, and non-standard components with 
volumes, such as rotor-stator structures. The former, to- 
gether with non-standard discrete losses, are modeled by 
Flowmaster [12-16], while the latter are modeled by 
Fluent, so that the flow conditions can be reflected in- 
stantly and more reliably. The interface between them 
transfers relative boundary conditions. Furthermore, Ab- 
aqus is applied to conduct a structural analysis and is 
coupled with Fluent, thereby allowing the interaction of 
the fluid field and structure to be taken into account. For 
the system shown in Figure 5, the position of the fasten- 
ing bolts influences the clearance gap of the labyrinth 
seals, as well as the flow and heat transfer performance. 
Therefore, the interaction of the structure and the flow 
system must be analyzed. 

Figure 6 shows the structured grid of the flow field of 
the rotor-stator structures, and Figure 7 shows the finite 
element model of the rotating disk. 


Fig. 6 Structured grid of the flow field of the rotor-stator 
structures 


Fig. 7 Finite element model of the rotating disk 


Figure 8 shows the modeling scheme of the flow paths 
of the secondary air system. The radial positions of each 
component are aligned with the nodes connected to the 
inlet and the outlet. This paper takes the starting process 
of the aero-engine as an example and analyzes the per- 
formance in unsteady state conditions. The unsteady re- 
sults are verified and compared with experimental data 
and simulation results in the steady state. 

The standard structural components, including the foil 
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nozzles, annulus passage, round orifices, labyrinth seals, 
and non-standard discrete losses were modeled in Flow- 
master. Two structures of the rotor-stator were modeled 
in Fluent, and the structural rotating disk was modeled in 
Abaqus. MpCCI was used to build the topological rela- 
tionship between different physical fields. Figure 9 
shows the scheme of the 1D-3D-Structure coupling. 
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Fig. 8 Scheme of the flow paths of the secondary air system 
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Simulation and Results 


Steady State 


The boundary conditions in a steady state are shown in 
Table 2. Because the accumulation or dissipation of mass 
and heat need not be considered, the sealing clearance of 
the labyrinths can be taken as a constant. In this way, the 
influence of the structure can be neglected. The heat 
transfer takes place primarily on the wall of the cavities. 
The heat flow through the disk’s outer boundary to the 
mainstream gas is given by existing data. Other walls can 
be taken as adiabatic. 


Table 2. Boundary conditions under steady state 


Q (r/min) Boundary Press. (Pa) Temp. (C) 
Node 1 200114 374.20 
Node 7 198113 374.20 
11440 
Node 11 190751 376.67 
Node 18 117017 377.25 
Node 1 275048 414.90 
Node 7 272298 414.90 
15400 
Node 11 258766 414.90 
Node 18 130033 420.40 
Node 1 432901 478.00 
Node 7 428572 478.00 
19800 
Node 11 406025 485.31 
Node 18 171318 487.01 
Node 1 557501 517.40 
Node 7 551926 517.40 
22000 
Node 11 523774 526.36 
Node 18 219152 528.44 


Table 3 Pressure change of some important components 
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Simulations in a steady state were carried out accord- 
ing to the working conditions listed in Table 2. The com- 
parison of the pressure distribution between the results 
and experimental data is shown in Figure 10. The com- 
parison of the pressure change of some of the main 
components is shown in Table 3, while the comparison of 
the volumetric flow and temperature distribution is 
shown in Table 4. 

These figures and tables show that the results of the 
simulations are highly consistent with the experimental 
data. The largest error in the pressure change, approx- 
imately 10.1%, appears at the position of labyrinth seal 2 
because its ratio in the whole system is very low and 
therefore this error is acceptable. The comparison shows 
that the calculation model can describe the flow charac- 
teristics of different components accurately. The temper- 
ature changes most significantly at the position of the 
rotor-stator structure, and the higher the engine’s rota- 
tional speed, the greater the temperature increases 
through the rotor-stator. In the simulation in a steady 
state, the error is due primarily to the 1D calculation 
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Fig. 10 Pressure distribution under steady state 


Q (r/min) Pressure change(Pa) Laby seal 1 Rotor-stator 1 Rotor-stator 2 Discrete loss 
Test 10658 -1326 -915 41559 
11440 Calculation 10456 -1244 -1007 41785 
Error (%) 1.89 6.18 9.14 0.54 
Test 19140 -2936 -1656 78039 
15400 Calculation 19234 -2814 -1701 78493 
Error (%) 0.49 4.16 2.65 0.58 
Test 33342 -6598 -3120 146323 
19800 Calculation 33771 -7003 -3301 146133 
Error (%) 1.27 5.78 2.45 0.13 
Test 43256 -9699 -4547 189856 
22000 Calculation 43315 -10137 -4671 189424 
Error (%) 0.14 4.32 2.65 0.23 


74 


Table 4 Volumetric flow rate and temperature distribution 


Vol. Flow (m*/s) 
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Temperature (°C) 


Q (r/min) Position 
Test Cal Error (%) Test Cal Error (%) 
Outlet of cavity 1 0.0275 0.0293 6.55 376.6 378.3 0.42 
Outlet of cavity 2 0.00885 0.00835 5.65 3712 380.1 0.76 
ed Outlet of laby seal 2 0.00536 0.00518 3.36 374.2 374.4 0.05 
Inlet of nozzles 0.0417 0.0428 2.64 374.2 - - 
Outlet of cavity 1 0.0408 0.0432 5.89 419.4 421.9 0.60 
Outlet of cavity 2 0.0122 0.0117 4.10 420.4 424.3 0.93 
ware Outlet of laby seal 2 0.007 0.00667 4.71 414.9 415.2 0.07 
Inlet of nozzles 0.06 0.06157 2.62 414.9 - - 
Outlet of cavity 1 0.0628 0.0655 4.30 485.3 489.4 0.84 
Outlet of cavity 2 0.0184 0.0172 6.52 487.0 492.5 1.13 
oe Outlet of laby seal 2 0.0102 0.0094 7.84 478.0 479.2 0.25 
Inlet of nozzles 0.0914 0.0921 0.77 478.0 - - 
Outlet of cavity 1 0.078 0.0813 4.23 526.4 533.2 0.01 
Outlet of cavity 2 0.0227 0.0212 6.61 528.4 535.7 0.01 
22000 
Outlet of laby seal 2 0.0127 0.0119 6.30 517.4 519.0 0.003 
Inlet of nozzles 0.1134 0.1144 0.88 517.4 - - 


model and the codes themselves. Some simplifications 
may also contribute to the errors, such as the fact that the 
heat transfer in some discrete losses and the heat transfer 
between the rotor disk and the axis were not considered. 


Unsteady State 


Using the startup of an engine as an example, this pa- 
per also investigated the unsteady state performance of 
the secondary air system by the 1D-3D-Structure coupl- 
ing method, allowing the characteristics of flow, heat 
transfer, structure, and their interactions during the un- 
steady process to be predicted more reliably and accu- 
rately. Because the convergence and the speed of con- 
vergence during unsteady calculation are influenced 
greatly by the initial field, especially for the CFD model, 
the interactive iteration of boundary parameters may 
make the convergence of the calculation much harder. 
Thus, before the unsteady calculation is carried out, it is 
necessary to have the initial field in a steady state. 


Working conditions 


The unsteady factors that must be considered include 
rotational speed change of the engine, total pressure and 
total temperature at the inlet of the whole system, total 
pressure or static pressure at the outlet of the whole sys- 
tem, and so on. Figure 11 shows the speed curve of the 
engine during a 60s startup process. 

Accordingly, the boundary conditions of the secondary 
air system are shown in Figures 12 and 13. 


Heat transfer in the system occurs primarily at the wall 
surface of the rotating disk. The heat flow between the 
outer disk surface and the mainstream gas is set accord- 
ing to the existing data. Other surfaces are taken as adia- 
batic. For the calculation in an unsteady state, the phe- 
nomena of mass and heat accumulation or dissipation are 
the main causes of unsteadiness in large systems. More- 
over, the influence of structure must be considered for 
secondary air systems. Because there are many sealing 
structures in the system (for example, the labyrinth seal), 


its sealing clearance is as low as a fraction of a millimeter. 


Thus, the deformation of relative components may influ- 
ence the sealing efficiency considerably. 
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Fig. 11 Speed curve of the engine during starting process 
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Fig. 13 Total temperature curve of boundary nodes 


Results 


It can be seen from Figure 5 that the rotational parts of 
labyrinth seal 1 and labyrinth seal 2 are fastened on the 
rotating disk with a group of screw bolts. Thus, the radial 
deformation of the disk will influence the sealing clear- 
ance directly, and then further influence the flow charac- 
teristics of the entire system. Figure 14 shows the varia- 
tions in the clearance of labyrinth seal 1 and labyrinth 
seal 2 during the startup of the engine. The deformation 
of the disk is due to the change in its centrifugal force 
and the consequent increase of thermal stress caused by 
the radial temperature difference. During this 60s process, 
the clearance decreased by 0.12mm, which accounts for 
20% of the initial sealing gap, and its influence on flow is 
obvious. 

Under the influence of the decrease of the clearance, 
the volumetric flowrate of the rotor-stator structures, la- 
byrinth seal 2, and the foil nozzles is shown over time in 
Figures 15-18, respectively. In these figures, the results 
of 1D-3D-Structure coupling, 1D-3D coupling that does 
not take into account the deformation of the structure, 
and the steady state are all included. 

As can be seen from these results, if the structure is 
not coupled, then the clearance of the labyrinth seals re- 
mains constant. The scale of the system itself is too small 
to accumulate or dissipate much mass and heat, and so 
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the results of 1D-3D coupling match those from the 
steady state well. In fact, the increase in engine speed and 
the radial temperature gradient decreases the sealing clea- 
rance of the labyrinths. As a result, the flow resistance 
increases, and the flowrate decreases. In this system, the 
flowrate of the branches that contain labyrinth seal 1 and 
labyrinth seal 2 decreases clearly by comparison to the 
flowrate in the steady state. Also, the decreasing ampli- 
tude increases gradually with the decrease in the sealing 
gap. At 60s, the flowrate of these two branches decreases 
by 14.9% and 23.8%, respectively, compared to the 
steady state, while the gross flow decreases by 15.3%. 
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Fig. 14 Sealing clearance of labyrinth seals 
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Fig. 15 Volumetric flowrate of rotor-stator 1 
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Fig. 16 Volumetric flowrate of rotor-stator 2 
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Fig.17 Volumetric flowrate of labyrinth seal 2 
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Fig. 19 Temperature rise of rotor-stator 1 


Under the same heat transfer conditions, the decrease in 
air flow is indicative of the decrease in dissipated heat. 
The radial temperature gradient inside the rotating disk 
may increase, as does the air temperature. As is shown in 
Figures 19 and 20, at the end of the 60s process of engine 
startup, the temperature of the air from rotor-stator 1 and 
rotor-stator 2 increases by 12.33% and 14.02%, respec- 
tively, compared to that in the steady state. This indicates 
that the cooling efficiency of the secondary air is lower 
than predicted. Therefore, it is necessary to take into ac- 
count the interaction of structure, heat, and flow. 
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Fig. 20 Temperature rise of rotor-stator 2 


In addition to the axial sealing change caused by the 
change in radial force, if the system contains a radial 
sealing structure, the axial force must also be investi- 
gated. As shown in Figure 5, the sealing structure be- 
tween nodes 15 and 16 adopts a radial sealing pattern, the 
flow capability of which may be influenced greatly by 
the axial force and clearance. However, because the 
sealing structure is assembled by the fixed constraint 
method, the total axial force caused by the pressure 
change in rotor-stator 2 is not sufficient to overcome the 
constraint. Thus, the influence of the axial force on the 
system is very limited. 

In summary, in the traditional simulation method of a 
secondary air system in a steady state, flow and heat 
transfer are calculated successively and the influence of 
structure is separated as well. This is acceptable for pre- 
dictions in the steady state. However, for the unsteady 
state, only the 1D-3D-Structure coupling method can 
integrate the flow, heat transfer, and structure calcula- 
tions to provide the researchers with more reliable data 
with which to test this theory. 


Conclusion 


This paper summarizes the basic theory and realiza- 
tion method of fluid solid coupling and 1D-3D coupling. 
Based on the interface code MpCCI, Flowmaster, Fluent, 
and Abaqus were integrated. Thus, the flow, heat transfer, 
and solid network of a secondary air system may be inte- 
grated by the 1D-3D-Structure coupling method. First, 
1D-3D coupling was applied to the simulation of the 
system in a steady state. The results of pressure were 
consistent with the test data. The greatest error in the 
pressure change (10.1%) appeared at labyrinth seal 2. 
Because this accounts for a very small percentage of the 
pressure change in the whole system, the error is accept- 
able. Furthermore, the largest errors in flowrate and tem- 
perature between simulation and the test were 6.6% and 
1.13%, respectively. 
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The results showed that the 1D and CFD models de- 
scribe the flow and heat transfer of the corresponding 
components accurately, while the 1D-3D-Structure 
coupling method is applied to the simulation of the sys- 
tem in the unsteady state. The results showed that the 
centrifugal force of the rotating disk increased with the 
engine speed, as did the radial temperature gradient in- 
side the disk and the successive thermal stress. Thus, the 
disk deformed in the radial direction. The clearance of 
the labyrinth seals decreased, the flow resistance in- 
creased, and the flowrate decreased. At 60s, the flowrates 
of the two branches that are connected with the labyrinth 
seals decreased by 14.9% and 23.8%, respectively, by 
comparison to the data in the steady state. The gross flo- 
wrate of the whole system decreased by 15.3%, as did the 
cooling effect. 

During the unsteady state process, the temperatures of 
rotor-stator 1 and rotor-stator 2 increased by 12.33% and 
14.02%, respectively, compared to the data in the steady 
state. Thus, it is clear that it is necessary to introduce the 
1D-3D-Structure coupling method into the unsteady state 
simulation of a secondary air system. Especially for 
small-scale aero-engines, the phenomena of mass and 
heat accumulation or dissipation are very weak. The un- 
steady state characteristics are reflected mainly in the 
interaction of flow, heat transfer, and structure. The 1D- 
3D-Structure coupling method can predict the unsteady 
state process more reliably and accurately. 


Acknowledgements 


This study was supported by funds from National nat- 
ural science foundation of China (Grant No.51176004). 
The support is gratefully acknowledged. 


References 


[1] Cross H., (1936), Analysis of flow in networks of con- 
duits or conductors, URBANA, University of Illinois 
Bulletin, Illinois. 

[2] Kutz K.J., Speer T.M., (1994), Simulation of the second- 
ary air system of aero engines, ASME Journal of Turbo- 
machinery, vol.116, pp.306-315. 

[3] Yannick Muller, (2009), Integrated fluid network-ther- 
momechanical approach for the coupled analysis of a jet 
engine, In ASME Turbo Expo 2009: Power for Land, Sea 
and Air, Orlando, USA, ASME GT2009-59104. 

[4] Schallhorn P. A., Hass N. E., (2004), Forward looking 
pressure regulator algorithm for improved modeling per- 
formance within the generalized fluid system simulation 


[5] 


[7] 


[8] 


[15] 


[16] 


Simulation of Unsteady State Performance of a Secondary Air System by the 1D-3D-Structure Coupled Method 77 


program, In AIAA/ASME/SAE/ASEE Joint Propulsion 
Conference and Exhibit, Fort Lauderdale, USA, AIAA 
2004-3667. 

Salvo R.D., Deaconu S., Majumdar A.K., (2006), Devel- 
opment and implementation of non-Newtonian Rheology 
into the Generalized Fluid System Simulation Program 
(GFSSP), In AIAA/ASME/SAE/ASEE Joint Propulsion 
Conference and Exhibit, USA, AIAA 
2006-4869. 

McAmis R.W., Miller J.T., Burdette R.R., (1996), Mod- 
eling fluid flow networks, In AJAA/ASME/SAE/ASEE 
Joint Propulsion Conference and Exhibit, Lake Buena 
Vista, USA, AIAA 1996-3120. 

Benson R S., (1982), The thermodynamics and gas dy- 


Sacramento, 


namics of internal combustion engines[M]. Oxford: Cla- 
rendon Press. 

Morel T, Flemming MF, (1990), Characterization of ma- 
nifold dynamics in the Chrysler 2.2 S.I. engine by mea- 
surements and simulation, In SAE, Detroit, USA, SAE 
Paper 900679. 

Kolke A, Walhorn E, Hubner B, , Dinkler D, (2004), 
Strongly coupled analysis of fluid-structure interaction 
with free surface flow, Proceedings in applied mathemat- 
ics and mechanics.4(1), pp. 338-339. 

D.S.MILLER, (1996), Internal flow systems-2nd ed[M]. 
Cranfield, BHR Group Limited. 

Qun Zhang, Toshiaki Hisada, (2001), Analysis of flu- 
id-structure interaction problems with structural buckling 
and large domain changes by ALE finite element method, 
Computer Methods in Applied Mechanics and Engineer- 
ing. 1(1): 231-236. 

W.F.McGreehan, M.J.Schotsch, (1988), Flow characteris- 
tics of long orifices with rotation and corner radiusing, 
Journal of Turbomachinery, 110(4), pp. 213-217. 

GJ.Yoo, R.M.C.So, B.C.Hwang, (1991), Calculation of 
developing turbulent flows in a rotating pipe, Journal of 
Turbomachinery, vol.113, pp. 34-41. 

Riccardo Da Soghe, Bruno Facchini, (2011), Analysis of 
gas turbine rotating cavities by a one-dimensional model: 
definition of new disk friction coefficient correlations set, 
Journal of Turbomachinery, vol.133, pp. 021020-1-8. 
Qinxue Tan, Jing Ren, Hongde Jiang, (2009), Predic- 
tion of flow features in rotating cavities with axial throu- 
ghflow by rans and les, In ASME Turbo Expo 2009: 
Power for Land, Sea and Air, Orlando, USA, ASME pp. 
GT2009-59428. 

Hay N., and Spencer A., (1992), Discharge coefficients of 
cooling holes with radiused and chamfered inlets. Journal 
of Turbomachinery, vol.114, pp. 701-706. 


